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Abstract 

The integrability of one dimensional quantum mechanical many-body 
problems with general contact interactions is extensively studied. It is 
shown that besides the pure (repulsive or attractive) 5-function inter- 
action there is another singular point interactions which gives rise to 
a new one-parameter family of integrable quantum mechanical many- 
body systems. The bound states and scattering matrices are calculated 
for both bosonic and fermionic statistics. 
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Quantum mechanical solvable models describing a particle moving in a local singular 
potential concentrated at one or a discrete number of points have been extensively dis- 
cussed in the literature, see e.g. [I], ||, [| and references therein. One dimensional problems 
with contact interactions at, say, the origin (x = 0) can be characterized by separated 
or nonseparated boundary conditions imposed on the (scalar) wave function (p at x = 0. 
The classification of one dimensional point interactions in terms of singular perturbations 
is given in j|. In the present paper we are interested in many-body problems with pair- 
wise interactions given by such singular potentials. The first model of this type with 
the pairwise interactions determined by (^-functions was suggested and investigated in |5[] . 
Intensive studies of this model applied to statistical mechanics (particles having boson 
or fermion statistics) are given in || |7| (these also leads to the well known Yang-Baxter 
equations). 

Nonseparated boundary conditions correspond to the cases where the perturbed op- 
erator is equal to the orthogonal sum of two self-adjoint operators in L 2 (— oo,0] and 
L,2[0, oo). The family of point interactions for the one dimensional Schrodinger operator 
— -rTf can be described by unitary 2x2 matrices via von Neumann formulas for self-adjoint 
extensions of symmetric operators, since the second derivative operator restricted to the 
domain C£°(R \ {0}) has deficiency indices (2,2). The boundary conditions describing 
the self-adjoint extensions have the following form 

(J?l-(:i)(J?l- 

where 

ad — bc—1, 9,a,b,c,d e M. (2) 

tp(x) is the scalar wave function of two spinless particles with relative coordinate x. (]T|) 
also describes two particles with spin s but without any spin coupling between the particles 
when they meet (i.e. for x — 0), in this case (p represents any one of the components 
of the wave function. The values 6 = b = 0, a = d = 1 in ([[J) correspond to the case 
of a positive (resp. negative) 5-function potential for c > (resp. c < 0). For general 



2 



a, b, c and d, the properties of the corresponding Hamiltonian systems have been studied 
in detail, see e.g. |4], §, [J. 

The separated boundary conditions are described by 

<p'(0+) = /iV(0+) , ¥>'(0_) = h-<p(0J), (3) 

where /i 1 * 1 G M U {oo}. /i + = oo or h~ = oo correspond to Dirichlet boundary conditions 
and h + = or far = correspond to Neumann boundary conditions. In this case it is 
impossible to express the perturbed operator as the orthogonal sum of two self-adjoint 
operators in L 2 (— oo, 0] and L 2 [0, oo). 

In the following we study the integrability of one dimensional systems of iV-identical 
particles with general contact interactions described by the boundary conditions ([!]) or 
(^D that are imposed on the relative coordinates of the particles. We first consider the 
case of two particles (N = 2) with coordinates x±, x 2 and momenta k%, k 2 respectively. 
Each particle has n-'spin' states designated by s± and s 2 , 1 < s« < n. For x\ ^ x 2 , these 
two particles are free. The wave functions if are symmetric (resp. antisymmetric) with 
respect to the interchange (xi,Si) <-> (x 2 ,s 2 ) for bosons (resp. fermions). In the region 
X\ < x 2 , from the Bethe ansatz the wave function is of the form, 

ip = a 12 e l(klXl+k2X2) + a 21 e l(k2Xl+klX2 \ (4) 

where an and a 2 \ are n 2 x 1 column matrices. In the region x\ > x 2 , 

tp = (P 12 a 12 )e i{klX2+k2Xl) + (P 12 a 2l )e i{k2X2+klXl \ (5) 

where according to the symmetry or antisymmetry conditions, P 12 = p 12 for bosons and 
P 12 = —p 12 for fermions, p 12 being the operator on the n 2 x 1 column that interchanges 
Si <-> s 2 . 

Let k\ 2 = (ki — k 2 )/2. In the center of mass coordinate X = (x\ + x 2 )/2 and the 

relative coordinate x — x 2 — x%, we get, by substituting Q) and (|) into the boundary 

conditions at x = 0, 

f a 12 + a 2i = e td aP 12 (a 12 + a 21 ) + ie l9 bk 12 P 12 (a 12 - a 21 ), , . 

1 ik 12 (a 21 - a 12 ) = e w cP 12 (a 12 + a 21 ) + ie ie dki 2 P 12 (a 12 - a 21 ) ^ ' 



for boundary condition ([[]), and 



ik 12 (a 21 - a 12 ) = h + (a 12 + a 21 ) , 
ik 12 P 12 (a 12 - a 21 ) = h_P 12 (a 12 + a 2 i] 



(7) 



for boundary condition @ respectively. 

Eliminating the term P 12 ai 2 from (||) we obtain the relation 

«2i = Y 2 12 a 12 , (8) 

where 

yl2 = 2ie ie k 12 P 12 + iknja - d) + (k 12 ) 2 b + c 
ik 12 (a + d) + (k 12 ) 2 b - c 

We remark that the system (|7]) is contradictory unless 

h+ = -h_ = heRU {oo}. (10) 

In this case it also leads to equation (§) with 

,12 ik 12 + h 



)-.r .'- . ■ (Mi 

>'l2 



21 ik„-h 



For > 3 and x\ < x 2 < ... < xn, the wave function is given by 

?P = a 12 ^ Ne ^x 1+ k. 2 x 2 +...+k N x N ) +Q , 2i ^ Ne i(k 2 x 1 +k 1 x 2 +...+k N x N ) + ( N \_ 2 ) other terms. (12) 

The columns a have x 1 dimensions. The wave functions in the other regions are 
determined from ( |T2"D by the requirement of symmetry (for bosons) or antisymmetry (for 
fermions). Along any plane Xi = sc»+i, i G 1,2, AT — 1, from similar considerations as 
above we have 

oti 1 i 2 ..j i i i+1 ...i N = Y%+ l .a h i 2 ... k+1 i i .a N -, (13) 

where 

= 2»e"fc, <t<+1 P" +1 + ik kk+1 {a -d) + (h lh+1 ) 2 b + c 

ik kk+1 (a + d) + (k hk+1 ) 2 b - c 1 1 

for nonseparated boundary condition and 

\rii+l _ ikkk+i + h / 1f -\ 
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for separated boundary condition. Here = (k^ — ki i+1 )/2 play the role of spectral 

parameters. p n+1 = p 1 ^ 1 f or bosons and P vl+1 = —p 1 ^ 1 for fermions, with p n+1 the 
operator on the n N x 1 column that interchanges Sj Sj+i- 

For consistency Y must satisfy the Yang-Baxter equation with spectral parameter 



, 1C , i.e., 



ym,m+lym+l,m+2ytn,m+l •\rm+l,m+2ym,m+lym+l,m+2 

ij kj ki ki kj ij ) 



or 



•ymr-yrs-ymr -yrs-yrar-Y^rs (16) 

ij kj k% ki kj ij V / 

if m, r, s are all unequal, and 

•ymr-ymr "ymr-ysq yS'/'^mr (IT) 

^3 3^ ' ^3 ^3 

if m, r, s, q are all unequal. 

The operators Y given by (14) satisfy the relation flT7|) for all 0,a,b,c,d. However 



the relations ( [TBI) are satisfied only when 9 = 0, a = d and 6 = 0, that is, according 
to the constraint (0), 9 = 0, a = d = ±1, b = 0, c arbitrary. The case a = <i = 1, 
9 = b = corresponds to the usual 5-function interactions, which has been investigated 
in || [7[]. The case a = d = — 1, 9 = b = 0, which we shall refer to as 'anti-5' interaction, 
is related to another singular interactions between any pair of particles (for a = d = — 1 
and 9 = b = c = see [|, Associated with the separated boundary condition, the 
operators Y given by (|I5|) satisfy both the relations ( |TTD and ([H]) for arbitrary h. 

We have thus found that with respect to iV-particle (either boson or fermion) problems, 
altogether there are three integrable one parameter families with contact interactions of 
type 5, anti-<5 and separated one, described respectively by one of the following conditions 
on the wave function along the plane x^ = Xj for any pair of particles with coordinates Xi 
and Xj, 

<p{0+) = +^(0_), y/(0+) = ap{0.) + y/(0_) , c e R ; (18) 
y?(0+) = -<p(0-), ^'(0+) = c^(0_) - y/(0_) , c G ffi ; (19) 
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y/(0 + ) = h<p{0+), <p'{0-) = -h<p{0-) , h E M U {00} . (20) 

The wave functions are given by (|1^) with the a's determined by ( |1"3D and initial condi- 
tions. The operators Y in (|I3]) are given respectively by 



v ii+l _ h i+1 )P U+1 + C _ . . 



i(k h -k h+1 )-c 



W <~ i(k h -k li+1 ) + c ' ^ 



■■ , -, — hi ,, ) + 2h , , 

i(k k -k li+1 )-2h • 1 ) 



and 



Nevertheless, from ([21]) and (|22|) we see that if we simultaneously change c — > — c and 
pn+i _p««+i^ these two formulas are interchanged. There is a sort of duality between 
bosons (resp. fermions) with 5-interaction of strength c and fermions (resp. bosons) 
with anti-5 interaction of strength — c. It can be checked that under the "kink type" 
gauge transformation U = Yii>j sgn(x, — Xj), the N-boson (resp. fermion) 5-type contact 
interaction goes over to the N-fermion (resp. boson) anti-5 interaction. Therefore these 
two situations are in fact unitarily equivalent under a gauge transformation IA that is 
non-smooth and does not factorize through one particle Hilbert spaces. 



The integrable system related to the case ( P3| ) is not unitarily equivalent to either the 
5 or anti-5 cases. In fact their spectra are different (see the bound states below). In the 
following we study further the one dimensional integrable iV-particle systems associated 
with (|2l 



When h < 0, there exist bound states. For N = 2, the space part of the orthogonal 
basis (labeled by ±) in the doubly degenerate bound state subspace has the form, in the 
relative coordinate x — x<i — x\, 

<il> 2> ± = (9(x)±9(-x))e h W. (24) 

The eigenvalue corresponding to the bound states fl24|) is — h 2 . By generalization we get 
the 2 Ar ( Ar_1 )/ 2 bound states for iV-particle system 

ipN, e _ = a e _ ]J(6(x k - Xl ) + e kl Q{x x - x k ))e h ^>i lx '~ Xjl , (25) 

k>l 



where a e is the spin wave function and e = {e kl : k > I}; e M = ±, labels the 2 7V ( 7V - 1 )/ 2 -fold 
degeneracy. 

It can be checked that ipN >e satisfies the boundary condition (p0|) at Xi = Xj for any 
i 7^ j G 1, N. The spin wave function a here satisfies P* J a = e^a for any i ^ j, that 
is, = e^-a for bosons and p^a = —e^a for fermions. ipN,e is of the form (|T2D in each 



region. For instance comparing with ( |12D in the region x\ < x 2 ... < xn, we get 

ki = ih(N — 1), k 2 = k\ — 2ih, k 3 = k 2 — 2ih, k^ = —k%. (26) 
The energy of the bound state ipN,t is 

h 2 

E = N(N 2 - 1) . (27) 

3 

The scattering matrix can readily be discussed. For real k\ < k 2 < .../cat, in each 
coordinate region such as x\ < x 2 < ...xn, the following term in (|12|) is an outgoing wave 

Vw = a 12 ... N e k ^ + - +kNXN . (28) 

An incoming wave with the same exponential as (^8|) is given by 

i> in = [P lw P 2 ( JV - 1 )...]a Ar(Ar _ 1) ... 1 e fc -^ + - + ^ 1 (29) 

in the region xn < £zv-i < ■■■ < x\. The scattering matrix is defined by ip out = Sipi n . 
From (Jl3|) we have 

&12...N — L r 21 r 31 ••• r 7Vl \ a 2...Nl — ••• 

_ rvl2v 2 3 y( A '- 1 ) iv l[yl2y23 v^ -2 )^ -1 )] [V 12 1™ , — Q'„, , 

— [ Y 21 Y 31 ■■■ r Nl \l Y 32 Y 42 — Y N2 \ " • [ Y N(N-l)\ a N(N-l)...l = O «JV(iV-l)...l, 

where Y™^, is given by (p3|). Therefore 

Cj _ gfpNlp(N-l)2 plN _ g/^pl2l |p23pl2l |p34p23pl2l j-p(JV-l)JV p!2l 



Defining 



we obtain 



X,, = Y^P'i (30) 



5— [X 2 iX 3 i...Xjv-l][-X'32X42...Xjv-2]-"[-^'jV'(JV-l)]- (31) 



The scattering matrix S is unitary and symmetric due to the time reversal invariance of 
the interactions. < s' 1 s 2 ---s' N \S\siS2...spf > stands for the S matrix element of the process 
from the state (fciSi, k 2 s 2 , k^s^) to the state (kis[, k 2 s' 2 , k N s' N ). The momenta (|26|) 
are imaginary for bound states. The scattering of clusters (bound states) can be discussed 
in a similar way as in 0. For instance for the scattering of a bound state of two particles 
(xi < x 2 ) on a bound state of three particles (x 3 < x 4 < x 5 ), the scattering matrix is 
S = [X32X42X52] [X31X41X51]. 

We have extensively investigated the integrability of one dimensional quantum me- 
chanical many-body problems with general contact interactions. Besides the repulsive 
or attractive 5 and anti-5 function interactions, there is another integrable one param- 
eter families associated with separated boundary conditions. From our calculations it 
is clear that these are all the integrable systems for one dimensional quantum identical 
many-particle models (of fermionic or bosonic statistics) with contact interactions. Here 
the possible contact coupling of the spins of two particles are not taken into account. A 
further study along this direction would possibly give rise to more interesting integrable 
quantum many-body systems. 
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